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^ I Motivated by recent 3-D calculations of the explosion of Type la supernova 

via a pure deflagration we calculate the observed spectra at 15-25 days past 
maximum light of a parameterized model which has a considerable fraction of 
unburnt C+0 in the central regions. Rather than attempting a self consistent 



3-D calculation, which is beyond the scope of current computer codes, we modify 
the composition structure of the 1-D deflagration model W7. In our exploratory 



O 

parameterized calculations, we find that a central concentration of C+O is not 
ruled out by observations for the epochs we study. We briefly examine whether 
nebular phase spectra could be incompatible with observations. 

X 

Subject headings: stars: atmospheres — supernovae: general 



1. Introduction 

Interest in Type la supernovae has been high over the last twenty years, since it was 
recognized that they closely meet the astronomer's criteria for a "standard candle". This 
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interest has grown with the realization that one can use the observed light curve to make a 
luminosity correction (similar to the Cepheid P-L relation, Phillips 1993; Riess et al. 1995, 
1996; Phillips et al. 1999). SNe la played an important role in current determinations of the 
Hubble constant with the Hubble Space Telescope (Sandage et al. 1996; Branch 1998) With 
the discovery of the "dark energy" (Riess et al. 1998; Perlmutter et al. 1999; Goldhaber et al. 
2001) enthusiasm for studies of Type la supernovae has grown. Heretofore, all calibrations 
of the luminosity of SNe la have been empirical in nature; relying largely on the sample 
obtained by the Calan-Tololo survey (Hamuy et al. 1996). While significant efforts are 
underway to improve the local sample of SNe la, theoretical understanding of the explosion 
mechanism and progenitor system are required in order to produce confidence in the validity 
of the empirically based results. 

In spite of the fact that the progenitor system for SNe la is still uncertain (Branch 
et al. 1995), there is little doubt that it requires the thermonuclear explosion of a C+O 
white dwarf. Given this fact, progress is being made in the understanding of the explosion 
mechanism itself and calculations have progressed from 1-D and 2-D, to full 3-D (Khokhlov 
1993, 2001; Hillebrandt & Niemeyer 2000; Reinecke et al. 1999, 2002c,a,b; Gamezo et al. 
2003) 

The calculation of the explosion mechanism is a formidable numerical problem. There is 
good agreement that the explosion begins as a deflagration with a laminar flame, which then 
develops into a fully turbulent flame via various instabilities, dominated by Rayleigh- Taylor 
instabilities. Even in 3-D, the flame cannot be fully resolved and a variety of subgrid models 
have been developed. Recent work has focused on calculations of the deflagration of C+O 
Chandrasekhar mass white dwarfs. While some calculations have had trouble obtaining the 
observed energy of the explosion, there seems to be general agreement that a fully resolved 
3-D calculation can produce significant kinetic energy and enough nickel to reproduce the 
general features of a Type la explosion. Recently, Gamezo et al. (2003) performed a series of 
3-D deflagration model calculations in which they obtained a total kinetic energy of about 
10 51 ergs. The key feature of the simulations is the highly convoluted turbulent flame sur- 
face that allows extensive interpenetration of burnt and unburnt materials. Similar results 
were found by Reinecke et al. (2002a,b). However, the strong interpenetration of burnt and 
unburnt material leaves a significant amount of unburnt C+O in the central regions of the 
ejecta. This result led Gamezo et al. (2003) to conclude that deflagration models would be 
in conflict with observed supernova spectra because there would be strong carbon or oxygen 
lines produced by the model, which are not seen in observed SNe la spectra. Based on this 
inference, Gamezo et al. (2003) rule out the deflagration model completely as a viable model 
for the explosion mechanism of SNe la. In drawing this conclusion Gamezo et al. (2003) rely 
on the work of Jeffery et al. (1992), Fisher et al. (1997), Mazzali (2001), and Kirshner et al. 
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(1993) which all use a highly parameterized direct spectral analysis technique for producing 
synthetic spectra. In particular, all these analyses make the Schuster-Schwarzschild approxi- 
mation and thus cannot address the formation of lines below the imposed photosphere. These 
analyses are extremely useful for line identifications in observed spectra and for determining 
the velocity interval at which certain ions must exist. However, these analyses cannot (and 
do not claim to) determine whether a given species must be absent. That is, from these 
analyses one can conclude that a certain ion must be present in a certain velocity interval, 
but one cannot conclude that any particular ion must be absent. 

In order to examine the true viability of a full 3-D hydrodynamical calculation, a full 
3-D radiative transfer calculation should be performed. However sophisticated 3-D radiative 
transfer codes do not yet exist that can tackle the problem of detailed radiative transfer in 
SNe la. In order to examine the viability of deflagration models we present the results from 
PHOENIX calculations using a parameterized version of the tuned 1-D deflagration model W7. 

2. Calculations 

The W7 model (Nomoto et al. 1984; Thielemann et al. 1986) is a good match to normal 
SNe la spectra (Harkness 1991a,b; Nugent et al. 1995a,b; Wheeler et al. 1995; Lentz et al. 
2001), especially the outer layers visible during the pre-maximum phase of the supernova. 
The outer layers of W7 are unburnt C+O enriched with solar abundance metals, expanding 
at velocities, v > 15000 km s -1 , with a total mass of ~ 0.07 M , followed by intermediate 
mass elements, with an interior composed mostly of 56 Ni but also including some other iron 
group elements. The model is a 1-D deflagration model, where the flame speed was adjusted 
to produce intermediate mass elements at high velocity as observed in SNe la spectra. 

We have altered the inner composition of the W7 to produce two alternate models with 
C+O mixed into the center. The models are prepared by removing a portion of the original 
mass in the W7 model and replacing it with an equal mass of C+O, with the C/O ratio 
= 1 by number, keeping the density structure the same. The C+O is used for the unburnt 
fuel mixed to the core found in 3-D deflagration models. In the first model, labeled "25% 
C+O core", we have replaced 25% of the inner core with C+O. The mixed material extends 
from the center to 3000 km s _1 , where the mixed fraction is then reduced linearly in velocity 
space to 0% at 4000 km s _1 . In the second model, designed to mock-up the Gamezo et al. 
(2003) model and labeled "50% C+O core", 50% of the mass is replaced with C+O from 
the center to 4000 km s -1 . The mixed fraction is reduced linearly in velocity space to 10% 
at 5000 km s^ 1 . The 10% C+O region extends to 15000 km s _1 where the W7 model 
becomes entirely unburnt material. Examining Figure 3 of Gamezo et al. (2003) it is diffiult 
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to determine the exact composition structure, since they only distinguish between unburnt 
(C+O) and "burnt" material, but the exact composition of the burnt material is unspecified 
because it is too computationally intensive to follow the reaction rates in detail. 

The spectra and 7-ray deposition are calculated using the general purpose radiative 
transfer code PHOENIX, version 11.9 (Hauschildt & Baron 1999). There is some change in 
the total gamma-ray deposition with the composition change, but the change is smaller than 
one might expect, since in this velocity range most of the material in W7 is neutron-rich iron 
peak elements and not 56 Ni. For the 25% C+O core model the change is 2%, and for the 50% 
C+O core model the change is 20% at both epochs. PHOENIX includes all of the effects of 
special relativity in the steady-state solution of the transport equation and energy balance. 
For these exploratory calculations we have treated all species in LTE. Studies by Baron 
et al. (1996) and Wheeler et al. (1998) have shown the LTE is a reasonable approximation in 
SNe la and since we are looking only qualitatively for lines that are clearly absent in observed 
spectra the assumption of LTE should not affect our conclusions. The use of PHOENIX with 
W7 models is described in detail in Lentz et al. (2001). 



3. Results 

Figure 1 displays the observed spectrum of SN 1994D on 1994 April 14, roughly 15 
days past maximum light (35 days past explosion, Patat et al. 1996) and the results of three 
calculations: 1) the original W7 model (denoted base), 2) the "25% C+O core" model, and 
3) the "50% C+O core" model. 

Figure 2 displays the composition structure of the three models at day 35. While C II 
lines become more prominent they do not necessarily rule out the existence of a substantial 
amount of C+O. Fig. 1 shows that as the C+O abundance is increased, the C II A4619 line 
becomes stronger, as does the C II A6580 line, and the C II A7234 line. In the UV there is 
no evidence for, e.g., C IV. In the "25 % C+O core" model, C+O mixing has not reached 
its peak until = 0.76 (the total continuum optical depth at 5000 A). 

To verify that we are indeed sampling the line forming region with our mixed models we 
calculated a model 44 days after explosion (roughly 24 days after maximum light) which is 
shown in Figure 3. In this case, r g ^ = 1 at 1100, 1450, 2000 km s _1 in the base model, the 
25% C+O model, and the 50% C+O core model respectively, since O serves as a source of 
free electrons. At this epoch the same features become too strong, although since our models 
are pure LTE and they use only the Kurucz (1993) line list which does not include forbidden 
lines. Therefore, lines of [Fe III] and [Co III] which are generally assumed to produce the 
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features near 5500 A and 6500 A (Bowers et al. 1997) are not included in our calculations. 
It is interesting to note that the velocity of the absorption minimum of the C II A6580 line 
in our model spectra is at about 6000 km s -1 , which is where C II is both prevalent, in our 
50% C+O core model, and is the top of the canonical line forming region Tg+j = 0.1. In 
the 25% C+O core model, the minimum moves somewhat redward to about 5400 km s _1 . 
However in this model there is no carbon between 4000 and 15000 km s _1 so the minimum 
is likely produced by a complicated blend of lines as is almost always the case in detailed 
models. The peak is almost certainly produced by C II A6580. 

The models of Gamezo et al. (2003) have both C+O mixing to low velocity, and also 
mixing of burnt material (possibly iron) to high velocity and we have not attempted to 
simulate that here. In addition the results of Gamezo et al. (2003) end at 1.9 s, before the 
explosion has become homologous, so it is difficult to translate their plot of compositions as 
a function of radial coordinate into a model in homologous expansion where both velocity 
and radius are equally good radial coordinates. 

4. Discussion 

Our results clearly demonstrate that one cannot simply examine the radial composition 
structure of a model and determine whether a hydrodynamical model is viable or not. On 
the other hand, we do claim to have shown that the model of Gamezo et al. (2003) is a 
viable model for SNe la. In particular, it is not possible to determine the exact composition 
structure from their paper; future work should attempt to make some educated guesses 
about the composition structure and calculate synthetic spectra directly from their (angle 
averaged) model. Sorokina & Blinnikov (2003) calculated multi-band light curves from the 
model of Reinecke et al. (2002b) and found that having burnt products at high velocity was 
somewhat helpful particularly in U and B. This is an indication that mixed models may be 
viable for SNe la, whether produced by deflagration or some other mechanism. Nevertheless, 
a full synthetic spectrum calculation is required before firm conclusions can be drawn about 
the viability or lack thereof of any particular hydrodynamical model. What we have shown 
is that unburnt C+O in the central regions are not necessarily ruled out at the epochs we 
study. 

One possible problem with the model of Gamezo et al. (2003) is that late time nebular 
spectra would show strong forbidden oxygen lines, which have not typically been observed in 
SNe la, although Minkowski (1939) did identify weak narrow lines in SN 1937C of [O I] A6300 
and A6394 beginning 184 d after maximum. The lines were absent by 275 d after maximum. 
Again, until detailed late-time synthetic spectra have been calculated, it is impossible to 
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draw a firm conclusion. It may be that since the material in the center of the Gamezo et al. 
(2003) model is dense and the radioactive 56 Ni is well mixed with the C+O that the oxygen 
stays ionized for a long time, but it does seem likely that [O I] lines should appear at some 
point. 

In addition, the models of Gamezo et al. (2003) mix burnt material into the outer layers 
as well as mixing unburnt material into the inner layers. We have not accounted for that 
effect in these simple exploratory calculations. However, we know that SNe la show a range 
of "temperatures", forming a "temperature sequence" that correlates with the light curve 
shape sequence (Nugent et al. 1995b, S. Bongard et al., in preparation). Exactly how the 
sequence is produced is not understood, although it is clear that models that produce more 
nickel will be both hotter and have higher opacity than those with less nickel (Khokhlov 
et al. 1993). "Hot" SNe la show lines of Fe III (Fisher et al. 1999; Hatano et al. 2002) 
and Si III (D. Branch et al., in preparation), but in parameterized models these lines are 
produced via extra thermalization (Nugent et al. 1995b) and not by direct excitation due to 
non-thermal electrons produced by 56 Ni decay. If very inhomogeneous mixing of burnt and 
unburnt material in the outer layers is a general feature of deflagration models (Khokhlov 
2001), then the observed uniformity of the Si II 6150A feature probably rules these models 
out (Thomas et al. 2002). 

From our simple parameterized results, we conclude that the viability or non-viability 
of any hydrodynamical model cannot simply be determined by the presence of a particular 
element in a certain velocity interval. The converse is not true, the absence of a particular 
element in a velocity interval where it is observed is fatal to any particular hydrodynamical 
model. It will be helpful when future hydrodynamical simulations can extend to the regime of 
homologous expansion and when the full compositions structure of the models are available. 
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Fig. 1 — The observed spectrum of SN 1994D obtained on 1994 April 5(Patat et al. 1996) is 
compared to synthetic spectra computed assuming this corresponds to 35 days past explosion 
using W7 compositions (base), W7 with C+O mixing only into the central regions (25% 
C+O core) and C+O in the core and intermediate velocity regions (50% C+O core). The 
line identifications are to guide the eye, the actual features are blends. 
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Fig. 2. — The unmodified (base) W7 compositions as a function of velocity are shown in 
panel (a), the core only mixed compositions (25% C+O core) are shown in panel (b), and 
the "50% C+O core" compositions are shown in panel (c). 
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Fig. 3. — The observed spectrum of SN 1994D obtained on 1994 April 14 (Patat et al. 1996) is 
compared to synthetic spectra computed assuming this corresponds to 44 days past explosion 
using W7 compositions (baseline), 25% C+O core, and 50% C+O core compositions. 



